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1. Introduction 

Precise lattice calculations of K ^ nn decays will provide quantitative insight into the origin 
of the A/ = 1/2 rule and direct CP violation in kaon decays. Previous calculations have relied 
on the quenched approximation or have attempted to use chiral perturbation theory to extrapolate 
from heavy quark masses down to physical masses, or both [|l|, ^, |3[ 0, ||, |6l, ^. The calculation 
presented here avoids both these sources of error by computing the K ^ nn amplitude directly 
using dynamical lattices with 2+ 1 flavours of domain wall fermions (DWF) at near physical pion 
mass. We use RBC/UKQCD 32^ x 64, Lg = 32 lattices which use the Dislocation Suppressing 
Determinant Ratio (DSDR) plus Iwasaki gauge action with inverse lattice spacing a^^ ^ lA GeV 
(j8 = 1.75) and domain wall height M5 = 1.8. We use ensembles generated with amf^ = 0.001, 
amf^'^ = 0.045, corresponding to a unitary pion mass of niji 180 MeV. 

2. Four- Quark Operators and The Effective Hamiltonian 

The weak interactions and the effects of heavier quarks can be included in the lattice QCD 
simulation by evaluating matrix elements of an effective Hamiltonian [|8|, ^]. In particular the 
conventions of [^] are used. We calculate matrix elements of four-quark operators between K and 
KK States. In this paper the amplitude Aj is calculated, which requires the evaluation of matrix 
elements of three operators. These operators are classified by how they transform under SU (3)l x 
SU {3)r: 2(27.1), 2(8.8) ^rid 2(8,8)mx- Pfogrcss in calculating Aq is described in [[l0|]. 

3. Boundary Conditions 

We wish to simulate the K ^ nn decay at physical kinematics, which requires the final state 
pions to have a non-zero momentum. This is achieved by imposing antiperiodic boundary condi- 
tions on the quark fields in one or more spatial directions. The allowed momenta of the quark are 
then given by pn = {k + Inn) /L, where L is the spatial extent of the lattice. 

We relate the physical \ Q%~Il'^2 \^^) matrix element to the unphysical matrix element 

(tt+tt+I Q%~^l'j2 1^^) using the Wigner-Eckart theorem. This simplifies the operators and allows 
us to use periodic boundary conditions on the up- and strange-quarks while using antiperiodic 
boundary conditions only on the down-quark, thus giving the two pions momentum while the kaon 
remains at rest. 

If antiperiodic boundary conditions are imposed on the d-quark in only the x direction with 
periodic boundary conditions in the y and z directions then we can have a two-pion ground state 
in which one pion has momentum Px = x ^^'^ other pion has momentum px = — The 
antiperiodic boundary conditions allow us to extract non-zero momentum pions without the need 
to fit to an excited state, which would have been necessary had we imposed periodic boundary 
conditions on all the quark fields. In principle we can impose antiperiodic boundary conditions on 
the d-quark in one, two, or all three spatial directions corresponding to individual ground state pion 
momenta of p = n/L, \/^n/L and \/^n/L respectively. 

4. Details of the Calculation 

The calculation was carried out on 62 configurations of dynamical 32^ x 64 lattices using 
DSDR-i-Iwasaki gauge action and domain wall fermions with Lg = 32, generated on BG/P machines 
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at Argonne National Laboratory. Further details of the ensemble generation are given in [11]. The 
inverse lattice spacing is a^^ = 1.365(22) GeV, the physical volume is (4.62 fm)^ and we set 
the light and strange valence quark masses to ami = 0.0001 and arris = 0.049 respectively. This 
corresponds to a pion mass of rrij^ = 141.9(2.3) MeV and a kaon mass of rrifc = 507.4(8.5) MeV. 

We combined propagators with periodic and antiperodic boundary conditions in the time di- 
rection in order to double the effective time extent of the lattice. The meson correlation functions 
contained propagators which were computed with a source at f = (corresponding to (P+A) be) 
and t = 64 (corresponding to (P-A) be). We also generated strange-quark propagators with sources 
at tK = 20, 24, 28, 32, 36, 40 and 44 in order to calculate K ^ Tin correlators with kaon sources at 
these times, while the two-pion sources remained at either ? = or ? = 64. Thus we could achieve 
time separations between the kaon and two pions of 20, 24, 28 and 32 in two different ways which 
doubled the statistics. These separations were chosen so that the signals from the kaon and two 
pions did not decay to noise before reaching the four-quark operator Q. 

For the kaon and pions with zero momentum we use propagators with Coulomb gauge-fixed 
wall sources. For the two pions with non-zero momentum we use the same type of propagators for 
the u quark but used propagators with antiperodic spatial boundary conditions for the d-quark with 
Coulomb gauge-fixed momentum wall sources of the "cosine" type 

■^p, cos (x) = COS ip^x) COS {pyj) COS {p^z) ■ (4. 1) 

We use the same cosine source for each d-quark, which introduces a cross term that couples to 
two-pion states with non-zero total momentum. For example, if we consider giving momentum in 
only the x direction the product of the sources of the two d-quarks is 



Sp, cos(xi)5p, cos(x2) = COS (^-^l) COS {^X2 



^ : e'r^ig'r-'^z +f>'f-'^if>-'r^2 +f^-'r-'^if>'r-^-2 +e-'r-^i--'r^2 



(4.2) 



We require the two pions to have individual momentum pi = jX and pi = — fx, but the first and last 



terms of equation (4.2) couple to two-pion states with total momentum 2£ and — 2£ respectively. 
We eliminate the unwanted terms in the two-pion correlator by using pure exponential momentum 
sinks which constrain the final state to have zero total momentum. In the K ^ nn correlator, the 
zero momentum kaon has a similar effect on the cosine sources of the two-pions. Had we used the 
more conventional momentum source 

^p(x)=e'P'' (4.3) 

we would have needed to perform two separate d-quark inversions with momentum +p for one 
and — p for the other. The cosine source eliminates one of these inversions. In practice we only 
compute d-quark propagators with antiperiodic boundary conditions in or 2 spatial directions, 
corresponding to pions with ground state momenta p = and p = ^/in/L. This choice is motivated 
by the expectation that for our choice of quark masses, p = V^n/L will correspond to on-shell 
kinematics. 

5. Analysis and Results 



We extract the K ^ nn matrix element ^ by fitting a constant to the left hand side of (5.1 ) 
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■t)Cn: 



(5.1) 



Cknn is the/: 



CK{tK-t)Cnn{t) ZrZ 

Tin 

Tin correlator with a kaon source at Ik, i labels the four-quark operator Qi which 
is inserted at time t, and Zk and Zj^k are calculated from the kaon and two-pion correlators respec- 
tively, whose sources are at f = 0. The left hand side of equation (5J ) is plotted in figure |l| for each 
of the three operators. The figure demonstrates that sufficiently far from the kaon and two-pion 
sources we are justified in fitting to a constant. The fit results for {ZxZt„^) are indicated on the 
plot. 
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(a) (27, 1) operator 
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(b) (8,8) operator 



(c) (8,8),nx operator 



Figure 1: — > %% quotient plots for p = ^/2n/L. The two pion source is at f = while the kaon source is 
at f = 24. The dashed line shows the error on the fit 



The finite volume matrix elements are related to the infinite volume amplitudes A, using the 
LeUouch-Liischer factor |12, 111]. In particular we have 



\AV4 



50 ^ 95 



(5.2) 



where the quantity in square brackets (denoted by LL in table ^ contains the effects of the LeUouch- 
Liischer factor beyond the free field normalization. Et^-^ is the energy of the two-pion state, 5 is 
the s-wave phase shift, v is a factor counting the free-field degenerate states, q-^ is a dimension- 
less quantity related to the individual pion momentum k-^ via q-^ = kjiL/ln and is a kinematic 



function defined in |12]. Once q^^ is known, 5 can be calculated using the Liischer quantisation 



condition [14|. 



mi = 5{kj,) + (^{q,,). 

EniL is found by fitting the quotient of correlators Cjtn/ {CnY ~ Ae^^' to extract AE 



(5.3) 



2Eji. We then get the two-pion energy by calculating AE + 2Eji, where in the case of /? = 0, 
Etc = fnji while for p = v^tt/L, Ej^ is found from a 2 parameter fit to the pion correlation function 
which also has p = \fl%lL. This method is preferred to directly extracting E-^j^ from the two-pion 
correlator because the quotient C-jinj iCjif' cancels some common fluctuations in the numerator and 
denominator and reduces the error 

The pion momentum k-,i in the two-pion state is determined from the two-pion energy using the 
dispersion relation Ej^ti = 2^^m\ + k'^. It differs from = 0, ^/1k/L due to interactions between 
the two pions. Results for Ejm, kji, qji and 5 are presented in table d(j)/dq can be calculated 
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analytically so the only unknown in equation (5.2) is 35 /dq. The results for the phase shift can be 
plotted against kj^ and compared with experiment 
good agreement with experiment. For p 



15|, |16p. This is done in figure [2(a)| and we see 
we make the approximation that 5 is linear with kj^ in 
we use the phenomenological curve 



order to calculate dS/dqji (see figure ^(b)[ ). For p 
[17] shown in figure |2(a)| to calculate the derivative of the phase shift at the corresponding value 
of qjt. The derivative of the phase shift is found to be a small factor in comparison with d<p/dqjz. 
Results for dip/ dq^c and d5/ dq^ are shown in table ^ 



Table 1: Two pion energy and s-wave phase shift 



p 


Enn (MeV) 


(MeV) 


qn 


5 (degrees) 





285.9(4.6) 
489.2(8.1) 


17.55(61) 
199.2(3.8) 


0.0655(21) 
0.743(11) 


-0.306(29) 
-10.4(3.3) 



Table 2: Contributions to Lellouch-Luscher factor 



p 


d(p/dqn 


d5/dq^ 


LL 





0.239(14) 
5.039(35) 


-0.0815(50) 
-0.2927(52) 


0.9636(22) 
0.933(11) 



This Calculation 
Hoogland et. al. (A) 
Hoogland et. al. (B) 
Losty et. al. 
■ Plienom. (Schenk) 
_l_ 



100 



200 
kp, (MeV) 



(a) Comparison of calculated phase shift with 
experimental results [|l^, [l6[p7|]. 
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(b) Zoom of 2(a) showing that for small k S is 
approximately linear. The scattering length is 
calculated using chiral perturbation theory mM- 



Figure 2: Plot of / = 2 two-pion s-wave phase shift against momentum kj^. The results from p 
p — \/2%/L are shown in red. 



and 



The amplitudes A; are related to the physical decay amplitude via 

Al = a-^ J K^y;, £Q(At)Z,;(At)A,-, (5.4) 

where C, are the Wilson Coefficients and Zij are the renormalization constants, calculated using 
non-perturbative renormalization (NPR). The factor 1^3/2 is needed to convert from the unphysical 
— )• 71+71+ amplitudes back to the physical — )• 71+ 7r" amplitudes. At present only Z(27 j) has 
been calculated; for the (8,8) and (8,8)mx operators (which mix under renormalization) we make 
the approximation Z,y = 0.9Z^5,j. A full calculation of Zjj for the (8,8) and (8,8)nix operators 
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on the lattice is cuirently under way [|T^]. Comparing Eju^ for p = ^/in/L with the kaon mass 
rriK = 507.4(8.5) MeV we see that the decay is nearly energy conserving, so we use the results 
from p = y/ln/L to compute A2. Results for Re(A2) and Im(A2) for the four different kaon source 
times are shown in table |3[ Our final result for Re(A2) and Im(A2) is an error weighted average 
(EWA) over the four kaon source times. 



Table 3: Final results for A2. The errors quoted are statistical errors only. 



tK 


Re(A2)(units of 10"** GeV) 


Im(A2)(units of 10"" GeV) 


20 


1.33(11) 


-8.11(52) 


24 


1.44(11) 


-8.77(60) 


28 


1.53(13) 


-8.58(58) 


32 


1.20(16) 


-9.01(75) 


EWA 


1.396(81) 


-8.46(45) 



6. Systematic Error 

The major sources of systematic error in the determination of A2 are scaling violations, finite 
volume effects, partial quenching, uncertainty in d5/dq, and the fact that the masses and momen- 
tum are slightly different from their physical values. Furthermore, the approximation made for the 
renormalization constants for the (8,8) operators introduces a large systematic error into Im(A2) 
which we will estimate as 20%. A2 is very sensitive to scaling violations because it is proportional 
to a^^. We estimate this systematic by calculating Re(A2) with a lattice spacing determined from 
/k, f^n, and ro respectively, and find a fluctuation of 8.5% among the three values. For finite vol- 
ume effects we estimate 7% for the systematic error using finite volume chiral perturbation theory 



for the K — KK matrix elements [20, 21]. One expects that for A/ = 3/2 decays partial quench- 
ing will introduce small errors and in [ |2^ the use of partial quenching has been shown to affect 
Re(A2) by about 2%. A value of d5/dq that is rather larger in magnitude is obtained just by putting 



a straight line through the two-pion phase shift data points from this calculation in figure 2(a)| ; this 



value yields a result for Re(A2) that differs by 2% which we use as our conservative estimate of 
this systematic. Finally, a K ^ nn calculation on 24^ quenched lattices was done for a variety of 
meson masses and two-pion energies [^], and shows that the deviations of these parameters from 
their physical values in the present calculation causes an 1.2% difference in Re(A2). Adding all 
errors in quadrature results in a preliminary estimate of 11% for the systematic error in Re (A 2) and 
23% for the systematic error in Im(A2). 

7. Conclusions 

We have presented preliminary results for the M = 2>/2 K ^ UK decay amplitude on 32^ lat- 
tices with 2 + 1 flavours of DWF and the Iwasaki-DSDR gauge action. We find = 141.9(2.3) MeV, 
rriK = 507.4(8.5) MeV and E^n = 489.2(8. 1) MeV. The main conti-ibution to Re(A2) is expected to 
be from the (27, 1) operator, and our result 1.396(081)stat(160)sys x 10^^ GeV can be compared to 
the experimental result of 1.5 x 10^^ GeV and is found to agree within error. This is the first time a 
calculation of this type has been achieved. Im(A2) is dominated by the operators in the (8, 8) rep- 
resentation, so we expect there to be a large systematic error on Im(A2) due to the approximation 
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made for Z^. This is reflected in our final answer Im(A2) = -8.46(45)stat(1.95)sys x lO'i^ GeV. 
This source of systematic error will be eliminated once the NPR calculation has been completed. 
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